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ABSTRACT 

We complement the PARSEC data base of stellar evolutionary tracks with new mod¬ 
els of massive stars, from the pre-main sequence phase to the central Carbon ignition. 
We consider a broad range of metallicities, 0.0001^ Z ^0.04 and initial masses up 
to Mini=35OM0. The main difference with respect to our previous models of massive 
stars is the adoption of a recent formalism accounting for the mass-loss enhancement 
when the ratio of the stellar to the Eddington luminosity, Pg, approaches unity. With 
this new formalism, the models are able to reproduce the Humphreys-Davidson limit 
observed in the Galactic and Large Magellanic Cloud colour-magnitude diagrams, 
without an ad hoc mass-loss enhancement. We also follow the predictions of recent 
wind models indicating that the metallicity dependence of the mass-loss rates becomes 
shallower when Pg approaches unity. We thus find that the more massive stars may 
suffer from substantial mass-loss even at low metallicity. We also predict that the 
Humphreys-Davidson limit should become brighter at decreasing metallicity. We sup¬ 
plement the evolutionary tracks with new tables of theoretical bolometric corrections, 
useful to compare tracks and isochrones with the observations. For this purpose, we 
homogenize existing stellar atmosphere libraries of hot and cool stars (PoWR, AT- 
LAS9 and Phoenix) and we add, where needed, new atmosphere models computed 
with WM-basic. The mass, age and metallicity grids are fully adequate to perform 
detailed investigations of the properties of very young stellar systems, both in local 
and distant galaxies. The new tracks supersede the previous old Padova models of 
massive stars. 
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1 INTRODUCTION 

Massive stars play an important role in the evolution of 
galaxies. They are the most important stellar sources of ion¬ 


izing and dissociation photons l 

Schaerer et al.|2011| 

Kimm 

& Cen|2014||Dale, Ercolano & Bonnell|2012 

|Gai et al.|2014| 

Hollenbach & Tielens[|1999| Yu, Wang & Li 

2015|). They in- 

ject a significant amount of kinetic energy through power- 


ful stellar winds (Mackey et al.|[2014|. They are among the 


main drivers of metal and dust enrichment in galaxies when 
they explode as core collapsed supernovae (SNe) (Sarangi 


I&: Cherchneff][^015| [Schneider, Ferrara fc Salvaterra||2004p . 

They are thus very important sources of feedback to the am¬ 


bient ISM (Dale, Ercolano & Bonnell 20131. Last but not 


least, most of our information on the ongoing star forma¬ 
tion rates across the Universe heavily relies on our detailed 


knowledge of their properties (Bruzual & Chariot 2003 Ken- 
|nicutt|1998| ). 

Because of their relevance for so many fields of astro¬ 
physics, massive stars have been the subject of many ob¬ 
servational and theoretical investigations that are impossi¬ 
ble to list here. Understanding their evolution is challenged 
by the complexity of several physical phenomena, as re¬ 


cently reviewed, for example, by Martins & Palacios (20131 


Padova (Bertelli et al. 20091, MESA ([Paxton et al. 


and STARE VOL (Decressinet al.|2009| evolutionary tracks 


who compare STERN (Brott et al. 20111, Geneva (Ek- 
[strom et al.[[2012[), FR ANEC ( [Chieffi fc Limongi[ 2013p , 


20111 


of massive stars. [Martins fc Palacios| ( [201^ find that, apart 
from the inclusion of rotation, the main differences among 
the models computed by different authors can be attributed 
to the different treatments of convection and mass-loss, with 
the former process being more important in the domain of 
the less massive stars. It is worth noting that the Padova 
tracks ([Bertelli et al.[[^009[) analyzed by [Martins fc Pala- 


(20131 are still those presented in Fagotto et al. (19941 


and Bressan et al. (19931. Since then, there have been many 
advances both in the basic input physics (opacities, equa¬ 
tion of state, nuclear reactions) and in the mass-loss theory 
(and observations), and these called for a substantial re¬ 
vision of the Padova code. But, while low and intermediate 
mass stars were systematically updated, with the last version 
being PARSEC (PAdova TRieste Stellar Evolution Code) 
models ( [Bressan et al. poT^ and refs, therein), massive stars 
have not been updated since then. 

In this paper we present the new evolutionary tracks 
of massive stars computed with PARSEC. Besides the ba¬ 


sic input physics, which is described elsewhere (Bressan 
[et al.[p012[ and refs, therein), the main novelties concern 
the recipes adopted for the mass-loss rates and the spec¬ 
tral energy distributions used to convert the tracks from the 
theoretical to the observational plane. Concerning mass-loss, 
there have been many efforts over the past years to deter¬ 
mine/predict the mass-loss rates of massive stars across dif¬ 
ferent spectral types and metallicities. It is now widely ac¬ 
cepted that hot massive blue supergiants (BSC) and Wolf- 
Rayet (WR) stars lose a prominent amount of their mass 
through line-driven stellar winds and that the mass-loss 


rates show a simple scaling law with the metallicity ( 

Gas- 

tor, Abbott & Klein||1975| Kudritzki & Puls 

20001 [Nugis & 

Lamers 2000 Vink, de Koter & Lamers 2000| 2001| [Mui- 

jres et al.|2012b 

Crowther|2007 Smithfl2014 

1 . In Luminous 


Blue Variable (LBV) stars, there is evidence that an im¬ 
portant mode of mass-loss is through eruptive mass-loss, 
that may contribute as much as or even more than the 
steady stellar wind (Smith 20091. In this phase th e mass- 
loss rate may easily reach a few of l0~‘^MQyr~^ (Lamers 


19891. This eruptive mechanism is still unknown but the 


observed rates could be explained by a super-Eddington 
wind ( Smith fc Owocld||2006 (, or by non-disruptive hydro- 
dynamic explosions (Barsukova et al. 20141. Interestingly, 


LBV stars are found near the so-called Humphreys-Davidson 
limit (Humphreys & Davidson 1979| ) which delimits the for¬ 
bidden region above which only very few stars are observed 
in the Hertzsprung-Russell (HR) diagram of the Galactic 
massive stars. In this respect, a remarkable result of recent 
investigations is that the mass-loss rates could be enhanced 
by a significant factor when the stars approach the Edding¬ 
ton luminosity ( Grafener fc Hamann|2008 [Vink et al.|2011 1, 
which is known to happen near the [Humphreys fc Davidson[ 
(19791 limit. For the later spectral types there are larger un¬ 


certainties both on the mechanisms and on the strength of 
the mass-loss rates. In red supergiants (RSG) one custom¬ 
arily adopts the observational parametrization by jde Jager,[ 
Nieuwenhuijzen & van der Hucht (19881, but the mass-loss 


rates in this phase are known to be uncertain by a large fac¬ 
tor ( Salasnich, Bressan fc Chiosi[[199'9 Meynet et al. 12015 1. 
As in less massive Asymptotic Giant Branch (AGB) stars, 
dust formation on the circumstellar envelopes could be one 
of the possible mechanisms responsible for this enhancement 
( van Loon et al.|2005 I. 

Concerning the atmosphere models used to predict the 
stellar magnitudes and colours, they are equally impor¬ 
tant for the interpretation of observed properties of mas¬ 
sive stars. For stars with negligible mass-loss, a compre¬ 
hensive stellar atmosphere library usually adopted is AT- 


LAS9 (Castelli & Kurucz 20041, consisting of plane paral 


lei models in local thermodynamic equilibrium. This library 
is particularly suitable for ~ A, F and G type stars. For 
cool giants, where the plane parallel and non-LTE (non¬ 
local thermodynamic equilibrium) approximations must be 
relaxed, comprehensive stellar atmosphere libraries are pro¬ 
vided by Phoenix ( Allard et al.|201^ and references therein) 
and MARCS ( Gustafsson et al.|2008 ( projects. For hot stars 
with high mass-loss rates a number of atmosphere mod¬ 
els have been released in the recent years like the WM- 
basic models ( [Pauldrach, Puls fc Kudritzki|19^ (, the CMF- 
GEN models ^Hillier fc Miller|1998 1 and the Potsdam Wolf- 
Rayet (PoWR) models (Grafener, Koesterke & Hamann 


2002 

IHamann & Grafener||2003 2004| Sander, Hamann & 

Todt 

2012 

[Hamann, Grafener & Liermann [2006 [Hainich 


et al. 12014 1. 


The paper is organized as following. In section 2, we 
describe our stellar evolution models for massive stars, with 
particular care to the description of the recipes adopted for 
the mass-loss rates. In this section we also compare the new 
models with our previous models and with one set among 


the most recent models found in literature, FRANEC (Chi 
effi fc Limongi|2013 1. In Section 3 we describe in detail the 
adopted atmosphere models and the procedure used to ob¬ 
tain as much as possible homogeneous sets of bolometric 
correction tables. Finally, in section 4, we discuss the re¬ 
sulting colour-magnitude diagrams predicted from the new 
evolutionary tracks and isochrones. 
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PARSEC evolutionary tracks of massive stars 3 


Our evolutionary tracks can be downloaded from 
http://people.sissa.it/~sbressan/parsec.html 
and the isochrones can be downloaded from 
http://stev.oapd.inaf.it/cgi-bin/cmd 


2 STELLAR EVOLUTIONARY TRACKS 

We compute new tracks of massive stars with the PARSEC 
code for a wide range of initial metallicities and for initial 
masses from Mini=14M0 to Mini=350 Mq. The evolution 
begins from the pre-main sequence phase and ends at central 
Carbon ignition. The new evolutionary tracks complement 
the already existing models of intermediate and low mass 
stars. 


2.1 Basic input physics 


The input physics used in PARSEC is already thoroughly 


described in Bressan et al. (20121, Bressan et al. (20131 


and Chen et al. (20141. Below we briefly summarize the 


main points with particular attention to those more rele¬ 
vant to the massive stars. The equation of state (EOS) is 
computed with the EreeEOS code of A. W. Irwiq^l Radia- 
tive opacities are from the OPAL project ( [Iglesias fc Rogers] 
1996|_ and, in the low-temperature regim e, from the /ESO- 
PURj project (Marigo & Aringer 20091. Conductive opac¬ 
ities are included following |Itoh et al.| ( |2008p . The nuclear 
reaction rates (p-p chains, CNO tri-cycle, Ne-Na and Mg-Al 
chains and the most important a-capture reactions includ¬ 
ing the a-n reactions) and the corresponding Q-values are 
taken from the recommended rates in the JINA reaclib data 
base ( Cyburt et al.||2010[). Electron screening factors for all 
reactions are from Dewitt, Graboske & Cooper (19731 and 


Graboske et al. (19731. Finally, electron neutrinos energy 


losses are computed following [Munakata, Kohyama fc Itoh 


(19851, Itoh & Kohyama (19831 and Haft, Raffelt & Weiss 


(19941. 


The metal abundances, O.OOOl^JZ^O.Od, follow the par¬ 
tition of [Grevesse fc Sauval| (|1998[ ) (GS98) with the recent 
revision by Caffau et al. (20111. With this revision, the cur¬ 


rent observed solar metallicity is Zq = 0.01524. The ini¬ 
tial Helium content at varying metallicity is determined by 
y = Yp -I- 1.78 X Z, with Yp = 0.2485 taken from |Komatsu| 
et al.| ( |201l| . 


2.2 Convection, overshooting and mixing 


The mixing length parameter (Bohm-Vitense 19581, de¬ 


rived from the solar model, is omlt = 1.74. We adopt the 
Schwarzschild criterion ( Schwarzschildj 1 S)58 1 to test the sta¬ 
bility of radiative zones against convection. In the presence 
of a gradient of chemical composition, an alternative crite¬ 


rion is that of Ledoux (Ledoux 19471. This condition may 


happen during the evolution of massive stars when the con¬ 
vective core grows in mass or when an intermediate radiative 
region of varying chemical composition becomes unstable to 
the convection. In this paper, we opt for the Schwarzschild 


^ http://freeeos.sourceforge.net. 

^ http://stev.oapd.inaf.it/aesopus. 


criterion because, on one side it has been shown that it is the 
more appropriate one to account for the effects of thermal 
dissipation ( Kato|[l966 1 and, on the other, the presence of 
a sizable overshooting region from the convective core sig¬ 
nificantly reduces the differences between models computed 
with the two alternative criteria ( Meynet fc Maeder||2000 1. 

Overshooting from the convective core is estimated 
within the framework of the mixing-length theory, allowing 
for the penetration of convective elements into the stable re¬ 
gions ( Bressan, Chiosi fc BerteIII||1981 h The adopted mean 
free path of convective elements across the border of the un¬ 
stable region, /c=Aci7p with A^ = 0.5, is calibrated on the 


colour-magnitude diagram of intermediate age clusters (Gl 


rardi, Rubele & Kerber|2009 1 as well as on individual stars 


(Kamath et al. 2010 Deheuvels et al. 2010 Torres et al. 


20141. We also account for overshooting at the base of the 


convective envelope, which is simply modelled by mixing 
the radiative region down to a distance of Ze=0.7Hp from 


the formal Schwarzschild border (Along! et al. 19911. We 


stress that the extent of the overshooting regions and the 
corresponding mixing efficiencies are still a matter of de¬ 
bate. Concerning core overshooting, a recent analysis of the 
period spacing of gravity modes in low mass Helium burning 
stars, suggests a quite sizable overshooting region, A^ = 1.0 
in the above formalism (Bossini et al. 2015, in preparation). 


Concerning envelope overshooting, work in progress (Tang 
|et al.| ( |2014[ ) and Rosenfield et al. in preparation) already in¬ 
dicates that using larger values of Ze (close to 2 Hp ) at the 
bottom of the convective envelope fits better the extended 
blue loops seen in metal-poor dwarf galaxies. Therefore, fu¬ 
ture releases of the PARSEC data base are likely to have 
these prescriptions revised. 

Finally, rotational mixing has not yet been introduced 
in PARSEC. 


2.3 Mass-loss rates 

At solar metallicity, the mass-loss phenomenon is known to 
dominate the evolution of stars for initial masses Mini^ 
30Mq. We account for this process following recent pre¬ 
scriptions found in the literature for the different spectral 
types. In the blue supergiant phase, for Ten Js 12000 K, we 
adopt the relations provided by |Vink, de Koter fc Lamei^ 
(2000 2001| . This formulation (Rvoi) shows an almost lin¬ 
ear overall dependence of the mass-loss rates on the metallic¬ 
ity, M oc (ZfMQjyr. In the supergiant phases with 
Teff < 12000 K we use the mass-loss rates provided by |de| 


Jager, Nieuwenhuijzen & van der Hucht (19881, Rdj, assum¬ 


ing the same dependence on the surface metallicity of Rvoi. 


For WR stars we use the Nugis & Earners (20001 formalism. 


RiVi. They also provide a dependence on the stellar metal¬ 
licity. The definition of the WR phases is provided below 


(see section 3.21. 

An aspect which is relevant for the more massive stars 
concerns the transition between the O-phase, the LBV/RSG 
phase and finally the WR phase and, most importantly, the 
dependence of this transition upon the metallicity of the 
stars. For example, in the old Padova models, e.g., |Bressan| 
et al. (19931, the transition to the super-wind phase corre¬ 


sponding to the LBV stars is artificially set at the stages 
when the models cross the Humphreys-Davidson instability 


limit (Humphreys & Davidson 19791 in the HR diagram. 
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Figure 1. Selected evolutionary tracks for massive stars with 
Z=0.02 (upper panel), Z=0.008 (middle panel) and Z=0.001 
(lower panel). The mass-loss rates are indicated with the colour 
bar. The thick black straight lines represent the Humphreys- 
Davidson limit ^Humphreys &; Davidson|1979[ | which delimits the 
forbidden region above which only very few stars are observed in 
the Hertzsprung-Russell (HR) diagram of the Galactic massive 
stars. The big solid and empty circles indicate the ZAMS and 
the end points of the tracks respectively. The triangles mark the 
beginning of WR phase. 


Figure 2. Upper and middle panels: comparison with previous re¬ 
leases of Padova evolutionary tracks (symbols, [Fagotto et al.|1994[ 
[Bressan et al.|1993[ l with Z = 0.02 (upper panel) and Z = 0.008 
(middle panel). Models of different masses are indicated with 
different colours and line styles/symbols. Note that the Helium 
abundances are Y = 0.28 for Z = 0.02 and Y = 0.25 for Z = 0.008 
in [Fagotto et al.j l [1994[ | ; [Bressan et al.j | |1993[ | . Lower panel: 
comparison with the FRANEC solar abundance (Z = 0.01345) 
models without rotation (symbols, jChieffi & Limongi|20T^ . The 
meaning of big circles is the same as in figure []J while the big 
triangles are used for the alternative models. The Humphreys- 
Davidson limit is also drawn as in figure^ 
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PARSEC evolutionary tracks of massive stars 5 


This is justified by the evidence that Galactic and Magel¬ 
lanic Clouds massive stars near this limit show mass-loss 
rates that may reach M ~10~^Mq/ut. However, while the 
Humphreys-Davidson limit is an observed property of the 
HR diagram of massive stars in near solar environments, 
it is used independently from the metallicity of the galaxy, 
in spite of the fact that the mass-loss rates themselves do 


of the effects of CNO abundances (Muijres et al. 2012a|. 


depend on the abundance of heavy elements ( 

Kudritzki & 

Puls 

2000||Puls, Springmann & Lennon|2000| |Mokiem et al. 

2007 

Smith||2014 

1 . This approximation becomes critical at 


very low metallicities. 

From the theoretical side, recent detailed studies of ra¬ 


diative wind models (Grafener & Hamann (2008 1 and Vink 
|et al.| ( [ 2011 [ )) show that the mass-loss rates are strongly en¬ 
hanced when the stars approach the electron scattering Ed¬ 
dington limit 


_ LKss _ ^ 

~ 4ncGM ^ 


( 1 ) 


Since for the most massive stars at solar metallicity this may 
happen near the Humphreys-Davidson limit, the above for¬ 
malisms could provide a modern description of the transition 


from 0-type through LBV/RSG-type to WR-types (Vink 
fc Grafener||2012 1 . We thus include in PARSEC the recent 
formulation of mass-loss rates by [Vink et al. (20111. The 
resulting HR diagram of a few selected evolutionary tracks 
for about solar metallicity is shown in the upper panel of 
figure For purposes of clearness, we do not plot the pre- 
main sequence phase. The colours along the tracks represent 
the strength of the mass-loss rates as indicated in the inset 
scale. In the figure, the thick black lines mark the so-called 
Humphreys-Davidson limit ( Humphreys fc Davidsoii]|1979 1 
which delimit the forbidden region above which only very 
few stars are observed in the HR diagram of the Galactic 
massive stars. As shown in the figure, the main sequence 
of the most massive stars extends up to the Humphreys- 
Davidson limit. Eventually, the stars may encompass this 
limit but the time spent in this region is very short because 
the mass-loss rate becomes so high that the stars rapidly lose 
their envelopes and turn into the hotter region of the HR 
diagram. We stress that this is not a result of an ad hoc as¬ 
sumption for the mass-loss rate, but a direct result of the ap¬ 
plication of the new adopted relations of the mass-loss rates. 
Near the Humphreys-Davidson limit, Fe rises close to 1 and. 


as described in Vink et al. (20111, when Fe is larger than 


0.7, the mass-loss dependence on Fe becomes high and, cor¬ 
respondingly, the mass-loss rates are significantly enhanced. 
With this formulation, the boosting of the mass-loss rate at 
the highest masses (MiniJS ISOM©) is effective already from 
the beginning of the main sequence and, because of the large 
mass-loss rates, they evolve almost downward vertically in 
the HR diagram. Interestingly, the luminosity of the tracks 
with the higher masses (Mini^ ISOM©) falls with time much 
more than those of the less massive ones. This is caused by 
the over-luminosity with respect to the main sequence mass 
- luminosity relation which, being larger at larger masses, 
results in larger values of Fe. Thus, the evolved brightest 
massive stars are not necessarily those with the largest ini¬ 
tial masses. 

The Fe dependence of the mass-loss rates of O-type 
supergiant stars has not yet been studied for different 
galactic environmental conditions, apart from the analysis 


A more thorough analysis in a broad metall icity range, 
10 ^ Zq Z ^ 2 Zq, has been performed by 


Grafener & 


Hamann (20081, but only for the case of WR stars. In partic¬ 


ular, Grafener & Hamann (20081 show that the dependence 


of mass-loss rates on the metallicity is also a strong function 
of Fe. While at low values of Fe the mass-loss rates obey 
the relation of M oc (.Z/Z©)° ®®M 0 /i/r, at increasing Fe the 
metallicity dependence decreases, and it disappears as Fe ap¬ 
proaches 1. In the absence of a more comprehensive analysis 
of the dependence of the mass-loss rates on the metallicity 
and Fe, and since there is a continuity between the models 
provided by [Vink et al. ( 2011| and those of WNL stars pro¬ 
vided by Grafener & Hamann| (2008 1 (see discussion in Vink 
erar]ponf )7 we assume in PARSEC that the scaling with 
the metallicity obeys the following relation 

M oc (Z/0.02)“ (2) 

with the coefficient a determined from a rough fit to the 


published relationships by Grafener & Hamann (2008 1 : 


Q; = 2.45-2.4*Fe 


(2/3 s: Fe < 1) 


(3) 


and with the supplementary condition 0 < a ^ 0.85. In 
summary, our algorithm for the mass-loss is the following. 
Besides the already specified mass-loss rate formulations 
(Rvoi, Rdj and Rjvi,) we compute also Rpe from the ta¬ 


bles provided by Vink et al. (2011 1 , but we scale the latter 
value with the metallicity, using equations it and §. Dur¬ 
ing the BSGs and LBVs phases the basic mass loss rate 
adopted is Rvoi. However, since, because of the effects of 
Fe, this can be encompassed by Rr^ and in order to secure 
a smooth transition, we adopt the maximum between Rvoi 
and Rpe- Toward the red supergiant phase we use Rdj, but 
again, since this is an empirical rate parameterized in such 
a way that it should hold on a broad region of the HR di¬ 
agram and it likely underestimates the mass loss rates of 
luminous yellow super-giants ( [Salasnich, Bressan fc Chiosi| 
19991, we compare it with Rr^ and take the maximum value. 


In the WR phases we consider the Nugis & Tamers (20001 
formulation. 

The HR diagram of a few selected evolutionary tracks 
for Z=0.008 is shown in the middle panel of figure 
We note here a significant decrease of the mass-loss rates, 
at a given mass. In particular while at Z=0.02 the mod¬ 
els of Mini~120M© and Mini~150M©, rapidly turn their 
main sequence evolution to higher effective temperatures, 
at Z=0.008 the mass-loss rates are not high enough to 
prevent the tracks from entering into the forbidden re¬ 
gion. Nevertheless, the tracks burn Hydrogen around the 
Humphreys-Davidson limit until, near central Helium igni¬ 
tion Fe becomes large and after performing a rapid excur¬ 
sion within the forbidden region, they turn into the blue 
part as WR stars. At even lower metallicities, the mass- 
loss rates decrease, unless the star is near the Edding¬ 
ton limit with Fe ~ 1 , and the location of the predicted 
Humphreys-Davidson limit shifts to higher luminosities. For 
example, at Z=0.001 in the lower panel of figure[^ the upper 
main-sequence widens significantly and the more massive 
stars evolve into the “forbidden” region even during the H- 
burning phase, because of their very large convective cores. 
They may also ignite and burn central Helium as “red” su¬ 
pergiant stars. 


© 2015 RAS, MNRAS 000,[iHT5l 
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2.4 Comparison with previous releases 


The current set of massive stars supersedes the old one 
adopted in several popular studies since 20 years (iBressan 


et al. 1993 Fagotto et al. 1994 hereafter BF+). With re¬ 
spect to BF+, the new models have a significantly finer 
mass spacing and extend up to a higher upper mass limit. 
The computed masses range from 14 to 20 Mq in steps of 
2 Mq , then up to 100 Mq in steps of 5 M© and finally 120, 
150, 200, 250, 300, 350 Mq . This allows for a quite better 
interpolation in ages and masses and a sampling of different 
initial mass function up to larger initial masses. Simple stel¬ 
lar populations can now be sampled with a higher accuracy 
than before and with more suitable mass and time-steps. 
The new tracks include also the pre-main sequence phase 
that begins when the central temperature of the protostar 
becomes larger than log(Tc/A') = 5.3. No mass accretion is 
accounted for during the pre-main sequence phase. 

To summarize the differences brought by the adoption 
of the updated physics input we compare the new tracks 
with those of the previous Padova release, in the upper and 
middle panels of figure for Z — 0.02 and Z = 0.008 re¬ 
spectively. Besides the presence of the pre-main sequence in 
the new tracks, which for purposes of clarity is not shown, a 
few general trends can be seen in the figure. At Z=0.02, the 
Zero Age Main Sequence (ZAMS) is similar. Note that the 
old tracks do not include the pre-main sequence phase. The 
main sequence termination is only slightly hotter in the new 
tracks, likely because of differences in the underlying opac¬ 
ities. The red supergiant phase is slightly cooler in the new 
tracks, but the differences are barely significant. We remind 
that in BF-I- the density inversion (arising in the inefficient 
region of the convective envelope) was inhibited in the com¬ 
putation of massive star tracks. In a later revision of the 


tracks (Girardi et al. 20001 it was inhibited at all masses. 


As discussed in Along! et al. (19931, a density inversion may 


develop in the external inefficient convection zones because 
of the requirement of the hydrostatic equilibrium in a re¬ 
gion with a large super-adiabatic real temperature gradient 
(assuming typical values of the mixing length parameter). 
This situation which should lead to a Rayleigh-Taylor in¬ 
stability or even a significant increase in the mass-loss rate, 
may give rise to numerical instabilities which preclude the 
computation of the track. To inhibit the density inversion 
one may use a mixing length parameter proportional to the 
density scale height {Hp) which, by rendering convection 
more efficient in the region where the density has a relative 
maximum, prevents a large super-adiabatic real tempera¬ 
ture gradient. Alternatively one may impose that the real 
temperature gradient is limited by Vt ^ , 

which is the choice made since BF-I-. This also simulates 
a more efficient convection and prevents the development 
of numerical instabilities. However since a larger convec¬ 
tive efficiency results in a hotter red giant track, and since 
BF-I- red giant tracks of intermediate and low mass stars 
compared well with the corresponding observed colours, we 
follow their method and inhibit density inversion only for 
the most massive stars. Thus the most massive red super¬ 
giants have effective temperatures that are slightly hotter 
than those obtained by allowing density inversion to occur. 
A detailed analysis of this effect is clearly needed, but likely 
it must also take into account effects of dusty circumstel- 



Figure 3. Upper panel: comparison of WM-basic spectra 
(smoothed) with different mass-loss rates. The effective tem¬ 
perature, gravity, radius and metallicity {Z = 0.02) are the 
same, but the black solid line is for M = 10“®Mq yr~^ and 
the red dotted line is for M = 10“^Mq yr~^. The ioniz¬ 
ing photon (with wavelength less than 912 nm) number ratio 
(Nionizing(M = 10-®Mq yr— l)/lVionizing(M = IO’^Mq yr-^) 
supposing the same luminosity) is 0.197. Lower panel: com¬ 
parison of WM-basic spectra (smoothed) with different metallici- 
ties. The effective temperature, gravity, radius and mass-loss rate 
(M = 1O“®M0 yr~^) are the same, but the black solid line is for 
Z = 0.02 and the red dotted line is for Z = 0.008. The ionizing 
photon number ratio (Aionizing)-^ = 0.02)/Nionizing(-Z = 0.008) 
supposing the same luminosity) is 0.994 (almost identical). All 
the spectra are normalized over the wavelength range from 0.84 
to 1.44 pm. 


lar envelopes which are known to affect the colours of RSG 
stars. 

The mass-loss rates adopted here are less efficient al¬ 
ready at Mini=30 Mq . At this mass, the star ignites Gar- 
bon as a red giant while in BF-I- the star is already moving 
toward the hot region of the HR diagram. This effect is 
more striking in the more massive stars, for those reaching 
the WR stages. Gomparing the two models of 60 Mq we see 
that while the H-burning and the central Helium ignition 
phases are pretty similar, the final WR phase is very dif¬ 
ferent. The mass-loss adopted for the WR stages in BF+ 
is significantly larger than that adopted here and the star 
ends with a luminosity which is about one order of magni¬ 
tude lower. 

At higher masses, Mini=100 and 120 Mq , there are sig¬ 
nificant differences already in the core H-burning. The mass- 
loss rates are initially higher in the BF-I- models, but when 
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the effects of Fe become important, the mass-loss rates be¬ 
come comparable and the tracks evolve along the similar 
path. There remain a large difference in the mass-loss rate 
adopted in the advanced WR stages. 

We stress again that in BF-I- models the mass-loss is 
arbitrarily enhanced when the stars approach or encompass 
the Humphreys-Davidson limit, while in the current models 
the mass-loss enhancement is a result of the photospheric 
conditions. In this respect, it is very interesting to note the 
similarity between the tracks of Mini=60 Mq that can be 
taken as representative of the fact that the mass-loss en¬ 
hancement is no more imposed as before, but it is instead a 
natural consequence of the photospheric conditions. 

The comparisons with the models of Z=0.008, typical 
of the Large Magellanic Cloud, confirm the previous pic¬ 
ture. The ZAMS is pretty similar, but now the effects of 
a lower metallicity become visible, since the normalization 
metallicity for the mass-loss rate has been assumed to be 
Z=0.02. In the old models, the metallicity dependence was 
slightly lower, M oc (^/0.02)°'®, and the mass-loss rates at 
Z=0.008 were about 40% larger than in the current mod¬ 
els, keeping other parameters fixed. Thus, the new models 
evolve at slightly higher luminosity until the effects of Fe 
become important. Again there is a striking similarity in 
the turnover of the effective temperature in the models that 
overcome the Humphreys-Davidson limit. As before, in the 
final WR stages the stars evolve at a significantly higher lu¬ 
minosity than that of the old models. We note that, while in 
the advanced WR stages the luminosities of the new models 
are significantly higher, the corresponding lifetime is more or 
less unchanged (~ 0.1 Myr for 100 M© models). This means 
that the new models should contribute far more to the hard 
ionizing photons than the old ones. This will be the subject 
of a forthcoming paper dedicated to the integrated proper¬ 
ties of young star clusters. 

For the comparison with models by other groups, we se¬ 
lect the evolutionary tracks without rotation computed with 
the FRANEC code ( |Chieffi fc Limongi|[M^ . The HR dia¬ 
gram is shown in the lower panel of figure The FRANEC 
models (dashed lines) are for a composition of ^ = 0.01345 
and Y = 0.265 while PARSEC models (solid lines) are 
for a slightly different chemical composition, Z = 0.014 
and Y = 0.273. They adopt a mixing-length parameter 
a = A/Hp — 2.3 and a core overshooting region of 0.2Hp 


(Chieffi & Limongi 2013 Martins & Palacios 20131. One 


of the major difference is the presence of big blue loops in 
the FRANEC 20,30 and 40 Mq models, while the big blue 
loops are absent in our models. In this respect, we note that 
the main sequence termination of our models is significantly 
cooler and more luminous than that of FRANEC models. 
This is particularly evident in the model of Mini =40 Mq 
and it is indicative of a larger central mixing during the 
main sequence phase of our models, which is known to dis¬ 
favour the occurrence of blue loops. The larger extension of 
the loops in the FRANEC models could also be due to their 
use of the Ledoux criterion for the definition of the borders 
of the intermediate unstable regions. It is well known that 
the Ledoux criterion inhibits the development of the inter¬ 
mediate convective regions within the Hydrogen profile, at 
the end of the H-burning phase, favouring a deeper pene¬ 
tration of the convective envelope when the stars move into 
the RSG region. This is an important issue deserving further 


investigation. Preliminary models computed adopting either 
the Schwarzschild or the Ledoux criterion in the Hydrogen 
chemical composition profile, confirm that the main cause of 
the inhibition of blue loops is the extent of core overshooting 
(or another kind of extended mixing) during main Hydrogen 
burning phase. With a slightly reduced core overshooting, 
also PARSEC models with the Schwarzschild criterion show 
extended blue loops (Tang et al. 2015 in preparation). 

We also see that the behavior of the tracks within the 
Humphreys-Davidson region is very similar. In both cases, 
the models abandon this forbidden region after losing a sig¬ 
nificant fraction of their mass. In both cases, the mass-loss 
must be enhanced, but the invoked mechanism are differ¬ 
ent: Chieffi & Limongi (20131 invoke the dust mass-loss en- 
hancement ([van Loon et al.|2005|), while we invoke the Fe 


mass-loss enhancement. 


3 STELLAR ATMOSPHERE MODELS 

We convert from effective temperature, gravity and lumi¬ 
nosity to magnitudes and colours using theoretical atmo¬ 
sphere models. Since one of our aims is to provide also 
isochrones in the observational plane, we need a full set of 
atmosphere models encompassing a wide range of param¬ 
eters, i.e. masses, evolutionary stages and metallicities. In 
the following, we describe the different atmosphere mod¬ 


els adopted. We begin with ATLAS models (Kurucz 1993 


Castelli & Kurucz 20041, which are the canonical models 


suitable for intermediate and low mass stars. We comple¬ 
ment this library with our new calculations of wind models 


for hot massive stars with the WM-basic code (Pauldrach, 


Puls & Kudritzki 19861, and new models of WR stars from 


the Potsdam group 1 

Grafener, Koesterke & Hamann||2002 

Hamann & Grafener 

20031 2004| Sander, Hamann & Todt 

2012 

Hamann, Grafener & Liermann 2006[ Hainich et al. 


20141. For the coolest stars we adopt the Phoenix models 


I Allard et al. 19971, already extensively discussed in our pre¬ 
vious paper (Chen et al.||2014 (. 


3.1 Intermediate mass and low mass stars 

The core of our library consists of the plane parallel ATLAS9 


models (Kurucz 19931. These LTE (Local thermal equilib¬ 
rium) models are well suited to describe the atmospheres of 
intermediate and low mass stars of spectral type between 
A and K. The most recent ATLAS9 models are those com¬ 
puted by [Castelli fc Kurucz| ( |2004p. ATLAS9 mod els are 
based on the solar abundances by Grevesse & Sauval (19981 


and make use of an improved set of molecular lines including 
TiO, H 2 O, HI-Hi and HI — H+. The model grids are com¬ 
puted for Teff from 3500 K to 50000 K, logg from 0.0 dex 
to 5.0 dex and [M/H]=-b0.5, -tO.2, 0.0, -0.5, -1.0, -1.5, -2.0, 
-2.5, -3.5, -4 and -5.5. We limit the use of ATLAS9 mod¬ 
els to the temperature range of 19000 K> Teff >6000 K. At 
higher temperature in general we need to consider models 
with mass-loss while at lower temperature the Phoenix mod¬ 
els are more appropriate. 


© 2015 RAS, MNRAS 000,[iHT5| 





































































Y. Chen et al. 


0 


M 

O 




Log(Teff) Log(Teff) 


Figure 4. Evolutionary tracks for massive stars with Z=0.02 (left panel) and Z=0.008 (right panel). Different colours represent different 
evolutionary stages: black dotted lines for stages precedent of WR phases, red solid lines for models using PoWR WNL-H50 {Z = 0.02) 
or WNL-H40 {Z = 0.008), blue dash-dotted lines for WN models, and brown dashed lines for WC models. The meaning of big circles is 
the same as in figure^ The overplotted vertical lines delimit the coverages of different atmosphere models as indicated in the plots. The 
Humphreys-Davidson limit is also drawn as in figure]^ 


3.2 Hot massive stars 

O, B stars. For temperatures typical of O and B stars 
(60000 K> Teff >19000 K) we have generated a new library 


of models nsing the public code WM-basic j Pauldrach, Puls 
fc Kudritzkijl986 b This allows us to take into account both 
the effects of extended winds and those of non-LTE, since 
both effects may significantly affect the emergent spectra 
of hot stars. We have generated new sets of stellar spec¬ 
tral libraries covering as much as possible the space of pa¬ 
rameters of our new evolutionary tracks of massive stars, 
i.e. effective temperature, gravity, metallicity and mass-loss 
rate. The temperature grid ranges from logT^g — 4.3 to 
5.0 in steps of AlogTeft = 0.025 dex. At each temperature, 
we adopt a step in gravity of Alog g — 0.5 dex, with upper 
and the lower boundaries that depend on the Teff. Indeed, 
the highest gravity is determined by the fact that the line 
driven radiation force cannot initiate the stellar wind, while 
the lowest value corresponds to stability problems when the 
models approach the Eddington limit. At each grid point 
we consider three different values of the mass-loss rates, 
M = 10“^, 10“®, and 10~^ MQ/yr, encompassing typical 
values for O, B stars. 

We note that, besides Teff, logg and stellar radius (ra¬ 
dius defined at a Rosseland optical depth of 2/3, R*, 2 / 3 ), 
there are additional input parameters in WM-basic, k, a 
and 5, referring most specifically to the structure of the wind 
as shown by Castor, Abbott & Klein ( 1975|). Thus, for any 


point in the (effective temperature, gravity) grid we suitably 
change the other stellar parameters (radius, luminosity and 
mass) as well as the wind parameters a and k in order to ob¬ 
tain three consistent atmosphere models with the required 
values of mass-loss rate, M = 10“^, 10~®, lO“®M 0 / 2 /r. 
These models are used to derive the broad band colours 


for M = 10~Mq yr~^, 98 for M — 10~^Mq yr~^ and 86 
for M = 10~^Mq yr~^. The metal partition adopted in the 
new models is the same as used in [Bressan et al.| ( |2012[ ) , 
to say |Caffau et al.] (2011|. As an example, we show in fig¬ 
ure the effects of changing the mass-loss rate at constant 
metallicity (upper panel) and those of changing the metal¬ 
licity at constant mass-loss rate (lower panel), for models 
with ZogTeff = 4.6 and logg = 3.5. In the upper panel, the 
black solid line refers to the model with M = 1 O“®M 0 yr~^, 
while the red dotted one is for M = 1 O“^M 0 yr~^. We 
can see that the Hell continuum is more absorbed in the 
model with higher mass-loss rate. We notice also that, as 
expected, some spectral features which appears in emission 
in the higher mass-loss rate model, turn into absorption in 
the model with a lower mass-loss rate. 

In the lower panel, the black solid line refers to the same 
model as in the upper panel, while, the red dotted one is now 
for a model at the same grid point but with Z = 0.008. This 
comparison shows that the effects of changing the metallicity 
on these hot spectra are less pronounced than those arising 
from the variation of the mass-loss rates. 

Wolf-Rayet stars. WR stars typically have wind den¬ 
sities one order of magnitude larger than those of massive 
O-type stars. Spectroscopically they are dominated by the 
presence of strong broad emission lines of Helium, Nitrogen, 
Carbon and Oxygen. They are subdivided into different sub- 
types, one with strong lines of Helium and Nitrogen (WN 
stars), another one with strong lines of Helium and Carbon 
(WC stars) and a third one with strong Oxygen lines (WO 
stars). 

To reproduce their spectra we make use of the most re¬ 
cent library of WR models computed by the Potsdam group 


(PoWR) (Grafener, Koesterke & Hamann 2002 Hamann 


and magnitudes for any given effective temperature, grav- 

& Grafener 2003| 2004| Sander, Hamann & Todt 

20121 

ity and metallicity, easily interpolating for the different val- 

Hamann, Grafener & Liermann||2006 

Hainich et al. 

2014 

1 . 


ues of the mass-loss rate actually used in the tracks. Our 
new library consists of about 300 models for each metal¬ 
licity. For example for a = 0.02, we computed 105 models 


They provide models with metallicities corresponding to 
those of WR stars in the Galaxy, the LMC and the Small 
Magellanic Cloud (SMC). The model grids are parame- 
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Figure 5. Colour-magnitude diagrams in HST/WFC3 broad bands for tracks of Mini=150MQ (upper four panels) and Mini=50MQ 
(lower four panels) with Z=0.02 (left panels) and Z=0.008 (right panels) respectively. The meaning of colours and line styles along the 
tracks is the same as in figure]^ The meaning of big circles is the same as in figure^ 
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Figure 6. Comparison of CMFGEN | |Hillier &: Miller|1998[ l WC 
and WN models with the same effective temperature Tgg = 
120000 K and metallicity Z=0.02. Also shown are the transmis¬ 
sion curves of three HST/WFC3 filters (F438W, F555W and 
F814W). Notice that the two strong emission lines HeII(4700) 
and CIV(5800) can greatly influence the magnitude/colours in 
related passbands. 


terized with Tefi,T =20 (the effective temperature at radius 
where the Rosseland optical depth is 20), and the trans¬ 
formed radius Rt (because models with the same set of such 
parameters show the same emergent spectrum as discussed 
in Schmutz, Leitherer & Gruenwald (1992[)). The trans- 


formed radius Rt is given by (Schmutz, Hamann & Wes- 
solowski|jl989 1: 


Rt = RY 


10 ‘^MQ'yr ^ 


2500 km 


M 


2/3 


(4) 


Since i?l =20 is nearly equivalent to the hydrostatic radius of 
the evolutionary tracks, this quantity combines stellar radius 
and mass-loss rate of the tracks with the terminal velocity. 
For the same reason Tefi,T =20 corresponds to the effective 
temperature of our hydrostatic models T*fj. As shown in 


Hamann & Grafener (20041, at large Rt (or thin winds). 


Teff,T=20 — Teff,T=2/3- (In the case of our WM-basic mod¬ 
els, Rt are always > 1.5, so there is no need to use Rt to 
match atmosphere models and theoretical tracks). As done 


Schmutz, Leitherer & Gruenwald (19921, we match the 


WR atmospheres to the evolutionary tracks by interpolat¬ 
ing T*fj and Rt. 

While the observational classification is quite well de¬ 
fined ( |Crowther|2007| |Growther, De Marco fc Barlow|1998p , 
it is more difficult to assign the WR subgroups along the evo¬ 
lutionary tracks. We use the following convention to identify 
them. WR stars are classified as WNL (late) when the sur¬ 
face Hydrogen mass fraction Xh is below a given threshold, 
XwNL. When Xh = 0, they are classified as WNE (early) 
if NC^C) < NC*N) and as WG if NC^C) ^ IV(^^N). 
We finally assign the WO subtype wh en the condition 
NC-^ C) + NC^O) > iV(^He) is fulfilled ( [Smith fc Maeder 
|1991| ). We then match our type assignments with the mod¬ 
els provided by PoWR, which adopt a fixed composition for 
any given subtype. For example, for the Galactic metallicity, 
in the PoWR library Xwnl=0.5 and the different subtypes 
have the following surface compositions: 

• WNL-H50 for Ah = A:wnl, A:He = 0.48, A:c = 1E-4, Xn = 
0.015, Ao = 0 and Xpe = 0.0014; 

• WNL for Xh = 0.2, Xhb = 0.78, Xc = lE-4, X^ = 0.015, 
Xo=0 and XFe = 0.0014; 


• W^NE for Xh = 0, Xne = 0.98, Xq = IE—4, X^ = 0.015, 
Xo = 0 and Xpe =0.0014; 

• WG for Xh = 0, Xhb = 0.55, Xc = 0.4, Xn = 0, Xq = 0.05 
and Xpe = 0.0016. 

For the LMG, the surface Hydrogen threshold is Xwnl= 0.4, 
while for the SMC they provide models with both Xwnl= 0.6 
(WNL-H60 models) and Xwnl= 0.4 (WNL-H40 models). 
Note that PoWR does not yet provide spectral models for 
WO stars, which, however, are quite rare objects especially 
at lower metallicities. 


3.3 Cool stars 


As reviewed by Allard et al. (19971, the atmospheres of 


cool stars are dominated by the formation of molecules and 
eventually at very low temperatures by dust condensation 
that can affect the spectral shape significantly. A suitable 
set of ID, static spherical atmosphere spectral models ac¬ 
counting for the above effects has been developed in the 
recent years and is continuously maintained by the Phoenix 
group ( Allard et al.|201^ I. Among the different suites of li¬ 
braries, we use the BT-Settl models which contain the most 
updated and complete stellar parameter grid and are well 


tested against observations (Allard et al. 20121. The full 


BT-Settl models are provided for 2600 K ^ Tgu < 50000 K, 
0.5 < log g < 6, and metallicities 0.000003 ^ Z < 0.04. We 
adopt such models for temperatures Teg 6000 K. These 
tables for cool stars have already been used in |Chen et al.| 
(20141, in the context of low and very low mass stars. 


3.4 Synthetic bolometric correction tables 


A concern for the libraries described above is that the dif¬ 
ferent sets {WM-basic, PoWR, ATLAS9 and BT-Settl) are 
calculated with different metallicities and different parti¬ 
tions of heavy elements. In order to obtain homogeneous 
spectral libraries, at least as far as the metallicity is con¬ 
cerned, we first calculate the global metallicity from the de¬ 
tailed abundance values provided by each group. We note, 
for example, that the BT-Settl models adopt a solar parti¬ 
tion ( jAsplund et al.||2009|) different from the one adopted 
in PARSEG (Caffau et al. 20111 and that, at low metal¬ 


licities, their partitions are a enhanced. On the contrary, 
for our WM-basic models, we use the same metallicities as 
in PARSEC. We then interpolate each set of spectra on the 
global metallicity grid defined by the PARSEC models. Since 
PoWR models are provided only for three typical metallic¬ 
ities, our Galaxy, the LMG and the SMC, we use Galaxy 
models for Z ^ 0.01, LMG models for 0.01 > Z > 0.006 and 
SMC models for Z ^ 0.006. Unfortunately this procedure 
may introduce some error because the partition of heavy el¬ 
ements is not the same in the different libraries. For example 
when we re-scale the [Caffau et al.j ( [2011[ ) and the [Grevesse[ 
& Sauval (19981 partitions to a metallicity of Z=0.02, we 
find a difference of ~17% in the Iron content and ~5% in 
the Oxygen content. However we verified that these varia¬ 
tions have negligible effects on the broad band colours of 
hot spectra. On the other hand this is the best one can do 
until a homogeneous set of model atmospheres (i.e. with the 
same partition of heavy elements) is computed on a range of 
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Figure 7. Isochrones of different ages, as indicated by the labels, are shown for Z=0.02 (left panel) and Z=0.008 (right panel). Note 
that at young ages the intermediate and low mass stars are still on the pre-main sequence (grey dotted lines). The Humphreys-Davidson 
limit is also drawn as in figure]^ 


parameters broad enough. The ranges of effective tempera¬ 
tures encompassed by different spectral libraries are marked 
by vertical lines in the HR diagrams shown in figure]^ for 
Z = 0.02 (left panel) and Z — 0.008 (right panel), respec¬ 
tively. Overplotted are selected stellar evolutionary tracks 
from Mini=20MQ to Mini=350 Mq. The combined atmo¬ 
sphere models, that share the same global metallicity of the 
corresponding evolutionary tracks as shown in the HR di¬ 
agram, provide a very good coverage in terms of effective 
temperatures and gravities. They allow for an optimal calcu¬ 
lation of bolometric correction (BCa) tables that are used to 
convert from theoretical to observational diagrams. The de¬ 
tails of this process are thoroughly described in |Girardi et al.| 
(20021 and Chen et al. (20141, and are not repeated here. 


While for the evolutionary tracks we only consider the pho- 
tospheric magnitudes, the corresponding isochrones also ac¬ 
count for the effect of possible circumstellar dusty envelopes 
following the dust calculation recipes described in |Marigo| 
et al. (20081 and Nanni et al. (2013 20141. 


4 DISCUSSIONS AND CONCLUSIONS 

We present new evolutionary tracks of massive stars for a 
broad range of metallicities, 0.0001 ^ Z ^ 0.04 and for 
initial masses up to Mini=35OM0. At super-solar metallic¬ 
ity, the models extend up to Mini=2OOM0 (Z=0.03) and 
Mini=15OM0 (Z=0.04), respectively. The new models com¬ 


plement the already published PARSEC data base ( Bres- 
san et aL]|2012 1 and supersede the old Padova evolutionary 
tracks of massive stars which are more than 20 years old. 
The stellar models are evolved from the pre-main sequence 
phase to the central Carbon ignition. The mass grid is very 
well sampled and it is fully adequate to perform detailed 
investigations of very young stellar systems both from the 
point of view of the resolved populations and from their in¬ 
tegrated properties. 

We revise the scheme adopted for including the mass- 
loss rate in the evolution of massive stars, by combining re¬ 
cent recipes found in the literature. Among the new recipes, 
particularly important is the enhancement of the mass- 
loss rate when the luminosity approaches the Eddington 


limit. We show that with the recent formulation by [Vink] 


et al. (20111, the models naturally reproduce the observed 


Humphreys-Davidson limit observed in the Galactic and 
LMC HR diagrams. In previous Padova models this limit 
was used as a threshold to enhance the mass-loss rates, inde¬ 
pendently from the metallicity. In this paper, the role of the 
metallicity is described by means of a simple recipe derived 


from the models presented in Grafener & Hamann (20081. 


The metallicity dependence of the mass-loss rate is now de¬ 
scribed by a power law with an exponent a which depends 
on Pe (equation 0 ) , the ratio between the stellar luminosity 
and the Eddington luminosity. When Pe approaches unity 
the metallicity dependence drops significantly, allowing for 
relatively high mass-loss rates also in stars of low metallic¬ 
ity. At lower values of Pe the usual exponent is recovered. 


a=0.85. While the models of Grafener & Hamann (2008 \ re¬ 


fer particularly to the WR phase, it has already been shown 
that there is a continuity between the new |Vink et al.| ( |2011[ ) 
formalism and the Grafener & Hamann (20081 results for 


WNL stars. The result of the new mass-loss formulation is 
that the Humphreys-Davidson limit shows a clear depen¬ 
dence with metallicity and it even disappears at very low 
metallicities. 

Compared to our previous models with the same metal¬ 
licity ( [Fagotto et al.|[l994[ [Bressan et al.||1993l ), the major 
difference concerns the WR phases. The new models evolve 
in the late WR stage, at logTen > 5.0, with luminosities 
which are significantly higher than those of the old Padova 
tracks (~ 0.9 dex in \og{L/L q))- These differences result 
from the different mass-loss recipes used in the advanced 
WR phases. Since the lifetimes of the models in this phase 
remain more or less the same, the new models should con¬ 
tribute much more to the hard ionizing photons. At lower 
metallicity, Z5C0.008, the new mass-loss formulation intro¬ 
duces significant differences even in the earlier phases. Our 
models compare well with the recent evolutionary tracks 
of non-rotating stars of similar metallicity, computed with 
FRANEC ( jChiefli fc Limongi|[M3t . 

Besides evolutionary tracks, we provide new tables of 
bolometric corrections that allow for the conversion from the 
theoretical HR to the observed colour-magnitude diagrams. 
For this purpose, we assemble existing atmosphere libraries. 
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Figure 8. Colour-magnitude diagrams in HST/WFC3 broad bands for isochrones of 1 Myr (upper four panels) and 2.5 Myr (lower four 
panels) with Z=0.02 (left panels) and Z=0.008 (right panels). The meaning of colours and line styles along the isochrones is the same as 
in figure 
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ATLAS9, PoWR, Phoenix and new atmosphere models cal- 
cnlated on purpose with the WM-basic code. We merge the 
different libraries with interpolation on a global metallic- 
ity scale, providing quite homogeneous tables of bolometric 
corrections, at several metallicities. 

An example of the colour-magnitude diagrams of the 
tracks of Mini=lOOM 0 and Mini=50MQ is shown in figurej^ 
for Z=0.02 (left panels) and Z=0.008 (right panels), respec¬ 
tively. In the figure, we highlight the different evolutionary 
phases with different colours and line styles. Black dotted 
lines indicate the evolution precedent of the WR stages, 
while the red solid lines are used to indicate the transition 
phase from the LBV phase to the late WR stars (WNL-H50 
for Z = 0.02 or WNL-H40 for Z = 0.008 in the PoWR nota¬ 
tion). The blue dash-dotted lines are used for the other WN 
stages and the brown dashed lines are for the final WC/WO 
stages. 

As can be seen in the figure, the optical colours show 
some evident jumps when the star type changes from WN 
to WC. This is likely to be caused by the appearance of 
different strong emission lines in different sub-types. This 
can be seen from figure where we compare the spectra of 
one WC star (black solid line) and one WN (red dashed 
line) star with the same effective temperature. We show 
the optical region sampled by the broadband filters F438W, 
F555W, and F814W of the HST/WFC3 system. To better 
clarify this point we have used in the comparison the high 
resolution CMFGEN spectra ( Hillier fc Miller|[l998 1 , which 
better show how the emission lines affect the spectra and 
the broadband colours. Both spectra show a strong emis¬ 
sion HeII(4700A) line, falling within the F555W passband 
and touching the border of the F438W passband. The dif¬ 
ferences between the WN and the WC spectral types are 
mainly on the contribution of the strong CIV(5800A) line 
within the F555W passband, in the latter type. Thus, the 
flux in the F555W passband is heavily enhanced in the case 
of the WC star with respect to that of the WN star, and 
correspondingly the sudden transition from the WN to the 
WC types is accompanied by a jump in the F438W-F555W 
and F555W-F814W colours. It is important to stress that 
the CIV(5800A) doublet is very sensitive to the adopted 
atmosphere parameters which somehow challenges its pre¬ 
dictability (private communication with Helge Todt, Wolf- 
Rainer Hamann and Gotz Grafener). Furthermore, we have 
to note that the transition between WN and WC spectral 
types is smoother in the tracks than in the spectra, because 
the latter are computed only at discrete values of the el¬ 
emental abundances, and this could enhance the effect, at 
least in terms of evolutionary speed. 

Using the evolutionary tracks we compute new 
isochrones of young stellar populations, with the same pro¬ 
cedure already described in [Bressan et al.| | |2012[ ). A few 
examples are shown in the theoretical HR diagram of fig¬ 
ure]^ from very young to very old ages and for Z=0.02 (left 
panel) and Z=0.008 (right panel), respectively. By means of 
the new bolometric correction tables, we convert theoreti¬ 
cal isochrones into observational magnitudes/colours, in the 
same way used for the evolutionary tracks. In figure we 
show the colour-magnitude diagrams of isochrones at very 
young ages, 1 Myr and 2.5 Myr, and for Z=0.02 and Z=0.008. 
The colour codings are the same as in figure but for the 
pre-main sequence which is drawn in light gray. 


A preliminary comparison of the new models with the 
colour-magnitude diagrams of star-forming regions in nearby 
low metallicity dwarf irregular galaxies, has already been 
performed in Tang et al. (20141. The new models of mas¬ 
sive stars will be used to compute the integrated proper¬ 
ties of young star forming regions (Obi et al. in prepa¬ 
ration). The full sets of evolutionary tracks can be down¬ 
loaded from http: //people. sissa. it/~sbressan/parsec. 
html The isochrones can be downloaded from http: //stev. 
oapd.inaf.it/cgi-bin/cmd 
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